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Abstract. We present some results of an on-going project aimed at studying a sample of 
Galactic H n regions ionized by a single massive star to test the predictions of modern generation 
stellar atmosphere codes in the H Lyman continuum. The observations collected for this study 
comprise the optical spectra of the corresponding ionizing stars, along with imaging and long-slit 
spatially resolved nebular observations. The analysis of the stellar spectra allows to obtain the 
stellar parameters of the ionizing star, while the nebular observations provide constraints on the 
nebular abundances and gas distribution. All this information is then used to construct tailored 
photoionization models of the H n regions. The reliability of the stellar ionizing fluxes is hence 
tested by comparing the photoionization model results with the observations in terms of the 
spatial variation across the nebula of an appropriate set of nebular line ratios. 
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1. Motivation 

In the last decades, a great effort has been devoted to the development of stellar 
atmosphere codes for massive stars. The advent of the new generation of nlte, line blan- 
keted model atmosphere codes, either plane-parallel (tlusty, Hubeny & Lanz, 1995), or 
spherically expanded (fastwind, Puis et al. 2005; CoSTAR, Schaerer & de Koter 1997; 
CMFGEN, Hillier & Miller 1998; WM-basic, Pauldrach et al. 2001) is already a fact. Each of 
them uses specific approximations in the calculation of the stellar atmosphere structure 
and spectral energy distribution (SED). 

This new generation of stellar atmosphere models produce quite different ionizing SEDs 
from the ones produces by the previous plane-parallel, hydrostatic models (either lte 
by Kurucz 1991, or nlte by e.g. Mihalas & Auer 1970). Some notes on this, and on the 
consequences on the ionization structure of H n regions, can be found in Gabler et al. 
(1989), Najarro et al. (1996), Sellmaier et al. (1996), Rubin et al. (1995), and Stasinska 
& Schaerer (1997). Although the new predictions seem to go in the right direction (viz. 
Giveon et al. 2002, Morisset et al. 2004) non-negligible differences can still be found 
between the various stellar codes (see e.g. Mokiem et al. 2004, Martins et al. 2005, Puis 
et al. 2005). 

Since the FUV range of the stellar flux cannot be observed directly, it is crucial to find 
indirect tests to constrain it. Ionized nebulae have many times been claimed as potential 
tools to check the validity of the emergent SED predicted by the stellar atmosphere 
models. Here we consider the possibility of using tailored models of simple, spatially 
resolved Galactic H II regions ionized by single massive stars to this aim. In this work we 
present preliminary results for two apparently spherical Galactic H n regions (M 43 and 
Sh 2-112) to illustrate the capabilities of the methodology we have developed. 
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2. Nebular candidates: the Galactic HII regions M 43 and Sh2-112 
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Figure 1. WFC images of the Galactic HII regions M43 (up) and Sh 2-112 (down), showing the 
location of the small apertures obtained from the long-slit used for the nebular spectroscopy. 
The ionizing stars of these nebulae are HD 37061 (B0.5 V) and BD +35 4216 (08 V), respectively. 
Note the proximity of M 43 to the Orion Nebula (M 42) in the upper right image. 



3. The observational dataset 



STELLAR SPECTROSCOPY NEBULAR IMAGING 



IDS@INT2.5 



H2400S (1^7500) 
X :(4000 - 4550] A 

range * J & 

X :[4500- 5050] A 

range 

H1B0OV (R=B500) 
/ w [6350 - 6750] A 

SNR - 200 - 300 



WFC@INT2.5m 

Narrow bond filters: 
H.i. HJi. [SII], [OIII] 

Spatial resolution: 
0.33 arcsec/pix 

FOV 33 x 33 orcmin 



" Nebular spectroscopy is obtained extractirig small size 
apertures from the long-slit 



LONG-SLIT NEBULAR 
SPECTROSCOPY 

ISIS@WHT4.2m 

R600& fS>. 2.25 A) 

^^[3880 - 5146] A 

R316R (5?. = 4.65 A) 

X : [5100 -7520] A 
X : [7300 - 9720] A 

Slit length*: 3.5 arcmin 
Slith width: 1.5 arcsec 



Table 1: Summary of the observational dataset characteristics 
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4. Detailed spectroscopic analysis of the ionizing stars: HD 37061 
and BD+45 3216 
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Figure 2. The stellar parameters of the ionizing stars were determined by means of the fit of the 
H , He i-ii observed profiles with FASTWIND synthetic profiles convolved with the instrumental 
and rotational profiles (see Tables 1, 2a-b). In the case of BD +35 4216 an extra-broadening (i.e. 
macroturbulence) was needed to properly fit the observed lines. 



HD 37061 (B0.5 V) 
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BD +45 3216 (08 V) 



m. = 9.07 + 0.01 
A, = 1.99 t 0.29 
4- 1.74! 0.17 Kpt 

=> M. = -*.t±0.5 
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7.2t 1.7 
logd/Le) = 4,95+ 0.2 1 
M/W Q = 24 + 12 
logQ(H°) = 48.4 



Tables 2a-b: Results from the spectroscopic analysis of the stars. Photometric data and dis- 
tances to the stars, used to estimate their stellar radii, luminosities and masses, are also indi- 
cated. 
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5. Comparing the ionizing SEDs predicted by the various stellar 
atmosphere codes 

As a first step, we compared the predictions — in terms of the shape of the ionizing 
SEDs — of the four modern stellar atmosphere codes CMFGEN, wm- basic, fastwind, 
and tlusty. To this aim we specifically calculated models with the various stellar codes, 
using the same stellar parameters (those indicated in Tables 2a-b), and the solar set of 
abundances from Grevesse & Sauval (1998). 
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Figure 3. Comparison of ionizing SEDs from CMFGEN, WM-basic, fastwind, and tlusty models 
with the stellar parameters indicated in Tables 2a-b. All SEDs have been re-mapped to the same 
A-grid for a better comparison. 



1.0 
f- 0.5 
o 0.0 
<i -0.5 



HD 37061 



■ S4 
i □ m 



□ a 



□ TL - CM 
A WH - CV 

a fw - cm 



q> =_o o o ° t lit 1 »* t t J ( ™_«\. IIS « 

woioz<zw(ji<zinoo<zioz 




inoioz<z0 



Figure 4. To quantify the differences showed by Figs. 3, we plot the quantity AlogQ(X +l ) 
for the ions indicated in those figures, using a scaie in which Alog Q(H °)=0. This type of 
diagram gives us an idea of how the shape of the SEDs compares in the H Lyman continuum. 
Alternatively, the number of ionizing photons for the various ionic species relative to the total 
number of H° ionizing photons can be compared (see its utility in Sect. 6). 



There are several ingredients which can affect the ionizing SED predictions: 

(a) The ionization structure of the stellar atmosphere 

(b) The blocking of emergent flux by the large amount of weak metal lines present 
in this spectral region (line blocking) and its effect over the ionization structure of the 
stellar atmosphere (line blanketing) 

(c) The presence of strong resonant lines (either in absorption or emission) 

(d) The presence of a velocity field in the stellar atmosphere 



All these effects must be taken into account when comparing the ionizing SED. Depending 
of the spectral region, the importance of each one of them can be more or less important. 
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6. The stellar - nebular connexion 

The state of ionization of a nebula depends, in first approximation, on the ionization 
parameter ({/), and the hardness of the ionizing radiation. Following Vilchez & Pagel 
(1988): 

n(X' +1 ) TT g^) tfcdv Q(X*) 

where n(X J+1 ) is the number density of the i-times ionized atom of element X, and 
H v is the stellar Eddinton flux. Therefore, a ratio n(X I+1 ) /n(X I ) is, to a first approach 
and for a given U, proportional to the number of photons able to ionize X 1 relative to 
that of the Lyman continuum photons. 

Therefore, we can use nebular line ratios involving lines from to successive ions of 
the same element (i.e. giving the ionization degree of the nebula) as constraints of the 
ionizing stellar SEDs at various wavelengths. These are indicated in Table 3. 



7. Some useful nebular line ratios 
N, diagnostic T e diagnostic 

1 1 6716 L J 6584 

Ionization degree 



[S III] 9069 Hel 667B [O III] 5007 Optical 
[5 II] 6731 Hp [O II] 3727+ 





[Ar III] 9.0 
[Ar II] 7,0 


[N HI] 57.2 
[N H] 121.7 


[S IV] 10.5 
[S III] 1B.7 


[Ne HI] 15.5 
[Ne II] 12.8 




Nebular 


abundances 




[Ne ETC] 3869 [S III] 9069 
Hp Hp 
Hel 667B [S H] 6731 


[0 III] 5007 
Hp 

[O II] 3727+ 


[Ar IH] 7135 

Hp 

[N II] 6584 



Table 3: List of some useful nebular line ratios 
indicating the nebular physical properties (T e , 
jV e , and ionization degree) and abundances. 
The color code is the same as the one used in 
Figs 5 and 6. 



Hp Hp Hp Hp 



8. Tailored photoionization models: probing the ionizing SED 
predictions 

We then constructed photoionization models for the nebulae, using the ionizing SEDs 
commented in Sec. 5. We used cloudy (v02.07.01, Ferland et al. 1998). Tables 4a-b 
summarize the characteristics of the CLOUDY models. The nebular abundances were ini- 
tially derived by means of the usual nebular techniques (i.e. direct method), and then 
fine-tuned to fit the observational constraints. 

Figs. 5 and 6 compare the photoionization model results (as a function of the projected 
distance to the central star) with the nebular observational data obtained from the various 
apertures. 
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M43 (CLOUDY models) 



SH2-112 (CLOUDY models) 



Stellar SEDs from CMFEEN. TLUSTV. WM-basic and FASTWIMD 
models with Teff =31 Kk logg=4.1 ]a g Q(l-l") = 47.1 
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5p*ier?col constant density models N Ebu lar a bunda nces { I og X/ H) 
N„ = 20 cm ! ( = N,[SII] } 
R - 1Q 1 * cm Dust free 



H*> -0.95. N = -4.50, O = -3.5b. 
Me ; -4.22, s = -5.16. Ar = -5.90 



Tables 4a-b: Input parameters used for the CLOUDY photoionization models. 

M43: The output from photoionization models using SEDs from the various stellar 
codes are very similar, except for the [S in] / [S n] line ratio (as expected from the com- 
parison of the quantity Q(S + ) /Q(H°) in Fig. 4). The observational constraints point 
towards predictions from the spherical models; unfortunately, the nebular spectrum of 
M43 is contaminated by an external nebular emission (as can be noticed from the be- 
havior of the [Olll]/H a , [Nein]/H Q , and [O in] /[On] line ratios in the outer part of the 
nebula), not allowing to achieve a firm conclusion. The behavior of these line ratios may 
be explained by assuming that the nearby 07 V star 9 1 Ori C (the ionizing source of M 42; 
see Fig. 1) is ionizing the nebular material in front of M 43. This contamination is clearly 
noticed when considering lines from O 2+ and Ne 2+ , since these lines are expected to be 
emitted by M43 itself only in a small region close to the ionizing star (due to the low 
temperature of HD 37061), but other lines may also be affected. 
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Figure 5. Comparison of the output from the photoionization models with the nebular 
spectroscopic observations for the case of M43. 
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Sh2-112: In this case, photoionization models show that nebular results are very 
different depending on the considered SED. It is interesting to see how results presented 
in Fig. 6 (concerning the nebular ionization degree constraints) and Fig. 4 are related. 
Note the importance of the IR line ratios in this case (except for the optical [O in] / [O II] 
line ratio, it is the only way to test the stellar SEDs for A ^ 430 A). 




Figure 6. As Fig. but for the case of Sh2-112. 
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9. A final note 

• We have shown that the methodology presented here (based on a combined stellar- 
nebular study of spatially resolved HII regions ionized by a single massive star) is very 
promising regarding the study of the reliability of the ionizing SEDs predicted by stellar 
atmosphere codes. 

• We have presented here results using simple dust free constant density spherical 
models representing the nebular gas, but we are also exploring the effect of considering 
dust, more complicated nebular gas distributions, and the possible escape of ionizing 
photons in the models. It is important to test the effect of these ingredients before 
proposing conclusions about the ionizing SED predictions. 
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